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The work by T. K. Ghosh and K. Machida [cond-mat/0510160 and Phys. Rev. A 73, 013613 
(2006)] on the sound velocity in a cyhndricaUy confined Fermi superfluid obeying a power-law 
equation of state is shown to make use of an improper projection of the sound wave equation. This 
inaccuracy fully accounts for the difference between their results and those previously reported by 
Capuzzi et al. [cond-mat/0509323 and Phys. Rev. A 73, 021603(R) (2006)]. In this Comment we 
show that both approaches lead exactly to the same result when the correct weight function is used 
in the projection. Plots of the correct behavior of the phonon and monopole-mode spectra in the 
BCS, unitary, and BEC limits are also shown. 
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In their recent study on sound propagation in an elon- 
gated Fermi superfluid in the BEC-BCS crossover, Ghosh 
and Machida 1] have reported a calculation of the sound 
velocity ui for the case of a power-law equation of state 
(EOS) as previously analyzed for bosons by Zaremba 
and for fermions by Capuzzi et al. Their result is 
found to be incorrect due to the use of an improper pro- 
jection of the sound wave equation. 

The eigenvalue equation for small-amplitude density 
modes (5nq(rj^)e"'^ is obtained by linearization of the hy- 
drodynamic equations around equilibrium and reads 

-V^- [noV^(a/./an|„^„fc,)] , (1) 

where no is the equilibrium density profile, co the fre- 
quency of the perturbation and q its wave vector along 
z. Equation reduces to Eq. (8) in Ref. I"] for a power- 
law EOS fi{n) ^ Cn<. 

To obtain the dispersion relation Loijq) for any value 
of q one must resort to the numerical solution of the 
eigenvalue equation ^ . A possible method to solve such 
an equation consists of expanding the eigenmodes Jn^in 
a complete set of basis functions (see e.g. Zaremba 
as 



(2) 



where a = (n^, m) labels the basis functions, with the 
radial number and m the number for the azimuthal an- 
gular momentum. By inserting this expression into Eq. 

and projecting the result onto an element of the ba- 
sis, a matrix representation of the eigenvalue equation is 
found, which is suitable for a numerical solution. This 
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procedure allows some freedom in the choice of the basis 
and of the projection, as long as these satisfy the bound- 
ary conditions. However, in order to obtain a standard 
eigenvalue equation of the form 



Av 



(3) 



one must choose a projection in which the basis is or- 
thogonal Q. For the basis functions adopted in Ref. 
cf. Eqs. (10) and (11), the projection must be performed 
with a weight function w(r) cx (1 — f^)"'''" where 
f = v±/R and 70 = I/7 — 1, R being the radius of the 
density profile. The orthogonality condition thus reads 



w{f ) 5n*^{v) SUa' (f ) d^r OC 6aa 



(4) 



Therefore, for Eq. (13) in Ref. [H to be correct the inte- 
grals defining the matrix Maa' must include the weight 
function w{r) and thus read 

X P^7"'l"'l)(2f2 _ 1) p(rro..\m\)^2p2 _ 

where P/j"'™"* are Jacobi polynomials. Equation ||SJ is 
what should have been used in Ref. instead of Eq. 
(14), where the weight function w{r) is missing and the 
constant A takes a different value since 6na{rj_) has not 
been normalized with w(r). A Fermi superfiuid with 7 = 
1 corresponds to a Bose-Einstein condensate (BEC) of 
molecules and has been previously analyzed by Zaremba 
for bosonic atoms. In this case w{r) = 1 and Eq. @ 
reduces to Eq. (14) in 0. 

To further analyze how the correct orthogonality con- 
dition affects the results, we have numerically solved the 
eigenvalue equation for sound propagation in a superfluid 
Fermi gas in the BCS, unitary, and BEC limits. Our re- 
sults for the two lowest frequency modes as functions of 
q are shown in Figs. ^andE] The lowest mode, shown 
in Fig. ^ is sound- like and has a phononic dispersion re- 
lation at long wavelengths. We observe that the slope of 



2 




FIG. 1: Dispersion relation uj{q) (in units of ui±) for the 
lowest-frequency (sound) mode as a function of qR with R the 
radius of the fermion density profile in the BCS limit. The 
dash-dotted, dashed, and solid lines correspond to fermions 
in the BCS, unitary, and EEC limits, respectively. The EEC 
limit corresponds to y — 0.25 in Eq. (26) in P]. 
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FIG. 2; The same as in Fig. Qfor the monopole mode. 



the dispersion relation, i.e. the sound velocity, is lower 
than that found in Ref. 1] and bends down as q increases. 
The first excited state, displayed in Fig.jSJ corresponds to 
a monopolar compressional mode that for g = is purely 
radial. Furthermore its frequency is known analytically 
at g = 1^ as Wo = \/l0/3 w_l in the BCS and unitary 
limits and as uiq — 2llJj_ in the BEC limit. Although the 
q = BEC limit of the monopole is correctly quoted in 
Eq. (8) of Ref. Q , it is not correctly depicted in the cor- 



responding Fig. 3, which is to be compared with Fig. [21 
in the present work. We also note that the effective mass 
associated to this mode is also different from that found 
by Ghosh and Machida. 

To obtain an analytical expression for the sound veloc- 
ity ui = du}{q)/dq\q=Q one can use Eq. (13) in Ref. [3 
for the lowest-frequency mode and expand it to first or- 
der in . This demonstrates that the off-diagonal terms 
in Alaa' do not enter the calculation of the sound veloc- 
ity, as pointed out by Zaremba for bosons. From the 
definition 101 one obtains Moo = 7/(7 + 1) and thus 



Ml = 



7 



274-2 



VF 



(6) 



with vf = a/2/i7M and jl the chemical potential, in 
agreement with our result previously obtained in Q]. 
This is the sound velocity that also Ghosh and Machida 
should have obtained if they had taken into account the 
weight function w(r). Their improper projection of the 
eigenvalue equation leads them to the pathological ex- 
pression ui = \/ {2 — 7)7/4tiF, which predicts no sound 
propagation for 7 > 2. The same problem affects the cal- 
culation of the effective mass mi, for the monopole mode 
(cf. Eq. (27) in Ref. pj), which once corrected is 



Mhui_ (2-^27)3/2(1 + 3^) 



2^ 



7(1-1-7-1-272) 



(7) 



An alternative and more direct procedure to evaluate 
the sound velocity is the one that we have outlined in 
0. The spatial dependence is eliminated from Eq. ^ 
above by integrating in the (x, y) plane. This yields the 
dispersion relation 



uj{q) = q J no dfi/dn\n=no Suq (fr± j j 5 



Tin d r± 



1/2 



(8) 

for a perturbation with J Suq d-r± ^ 0. Hence, the 
calculation of the sound velocity requires only the ex- 
pression of Sriq calculated at g = 0. By using Suq^o — 
{d^/dn\n=no)~^ we obtain 



Ul 



(^j^ Jnod'^r±j^ J {dfi/dn\n^no) ^d^r± 



1/2 



(9) 

This expression provides the exact velocity of sound 
propagation in cylindrically confined hydrodynamic gases 
with EOS /i(n), and for /L((n) oc n"^ leads to Eq. ©. 
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